Multiphoton microwave photoresistance in a high-mobility two-dimensional electron gas 
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We report on experimental and theoretical studies of microwave-induced resistance oscillations in a two- 
dimensional electron gas over a wide range of microwave intensities. We observe a distinct crossover from 
linear to sublinear power dependence of the oscillation amplitude and a concomitant narrowing of the oscillation 
extrema. To explain our observations we propose a theory based on the quantum kinetic equation at arbitrary 
microwave power. Taken together, these findings demonstrate a crucial role of multiphoton processes at elevated 
microwave intensities. 



Microwave-induced resistance oscillations (MIRO) ap- 
pearing in very high Landau levels of high-mobility two- 
dimensional electron gas (2DEG) have been discovered more 
than a decade ago.''^ Since then, a variety of other remark- 
able low-field transport phenomena have been identified and 
actively studied. Among these are phonon-induced^ and 
Hall field-induced^ " resistance oscillations, several classes 
of "combined" oscillationsr^"— zero-resistance states, ''^"^'^ 
several distinct zero-differential resistance states ,^^''"^'"^"" — 
and, most recently, a very sharp and strong photoresis- 
tivity peak near the second harmonic of the cyclotron 
resonance. ■'^^^ Microwave-induced conductance oscillations 
and zero-conductance states have also been realized for elec- 
trons on liquid helium. ^^-^^ 

Despite a huge body of experimental^^^^ and 
theoretical^^! work on microwave photoresistance, there 
are still many unsolved puzzles, such as activated temper- 
ature dependence at the oscillation minima, the role 
of microwave polarization,— and the effect of an in-plane 
magnetic field. Another outstanding issue, which is the 
main subject of this Rapid Communication, is the dependence 
on microwave intensity which to date was not investigated in 
detail. 

MIRO can originate from either the displacement 
mechanismr^i^^j^2i^2i^ stepping from the modification 
of impurity scattering by microwaves, or the inelastic 
mechanism ,^^i^^'^^'^^~ — owing to the microwave-induced 
non-equilibrium distribution of electrons. In the most 
widely-used approximation of weak microwave intensi- 
ties, both mechanisms lead to a simple expression for the 
photoresistivity, 

^ ~ -AX^neAiV) sin 27re, (1) 

PD 

where p/j is the Drude resistivity, e — bjjuj^, uj = 27r/ 
is the microwave frequency, uj^ — eB/m* is the cyclotron 
frequency of an electron with an effective mass to*, A = 
exp(— tt/cJcTq) is the Dingle factor, tq is the quantum life- 
time, A{V) = rfP is the reduced amplitude, rj is mechanism- 
dependent scattering parameter, and V is the dimensionless 
microwave power which, for active circular polarization, is 
given hy^ 

Vie) = - = F (2) 



Here, /S^, = (uiTem) \ t^™ = riee'^ /2^/e^eom*c, 2^/e^ = 
y/e + 1, e = 12.8 is the dielectric constant of GaAs, vp is the 
Fermi velocity, and £ac is the external microwave field. Equa- 
tion ^ thus predicts linear power dependence of the oscilla- 
tion amplitude and a power-independent phase, 5 = n — e^ = 
1/4, where n is an integer and e+ = e at the closest maximum. 

Experimentally, however, all studies, ^'^i*^^ except 
Ref. 20, found strongly sublinear power dependencei^^ These 
reports indicate that Eq. ([T), while used widely to analyze and 
interpret the experimental data, is not valid under typical ex- 
perimental conditions. At the same time, theoretical studies 
of the high intensity regime have only considered the limit 
of smooth disorder ''^•''^•''^'^^ However, it is well established 
both theoretically and experimentally that in a realistic high- 
mobility 2DEG sharp disorder plays a crucial role in many of 
the above mentioned non-equilibrium phenomena,—^ includ- 
ing MIRO,^ 

In this Rapid Communication we report on experimental 
and theoretical studies of microwave photoresistance in a very 
high-mobility 2DEG performed over a three orders of magni- 
tude range of microwave intensities. We observe a distinct 
crossover from linear to sublinear power dependence of the 
oscillation amplitude. This crossover is accompanied by nar- 
rowing of the oscillation extrema, which is reflected as a de- 
creasing phase. To explain our observations we propose a 
theory based on the quantum kinetic equation. Applicable 
at arbitrary microwave power, our model takes into account 
both the displacement and the inelastic contributions within a 
model of mixed disorder. Taken together, the findings demon- 
strate an important role of multiphoton processes at elevated 
microwave power and should be useful in interpreting prior 
experiments ^Vi^s- io 

Our Hall bar sample was fabricated from a symmetrically 
doped GaAs/AlGaAs quantum well with density rig — 2.9 x 
10^^ cm~^ and mobility ^ — 2.4 x 10'' cm^A'^s. The mi- 
crowave radiation of frequency / = 33 GHz was delivered to 
the sample using a WR-28 waveguide. Resistance was mea- 
sured at T ~ 1.5 K using a standard low frequency lock-in 
technique. Microwave intensity was varied by a calibrated 
attenuator Maximum power at the end of the waveguide is 
estimated as Pq ~ 1 mW. 

In Fig. [T] we present microwave magnetoresistivity Puj{B) 
measured at attenuations from 30 dB to 15 dB [Fig.fTfa)] and 
from 15 to dB [Fig.[nb)], in a step of 2.5 dB. At low inten- 
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FIG. 1: (Color online) Normalized microwave magnetoresistivity 
Pui{B) measured at attenuations (a) from 30 dB to 15 dB and (b) 
from 15 to dB, in a step of 2.5 dB. Integers mark MIRO order. 

sities [Fig.[Tfa)] the MIRO amplitude grows roughly linearly 
(5 dB is equivalent to a factor of about 3) and the positions 
of the maxima/minima remain constant (cf. I), in agreement 
with Eq. ([Til. At higher intensities [Fig.lTJb)] the amplitude 
grows considerably slower and eventually starts to decrease, 
likely as a result of heating.^"* At the same time, the extrema 
become sharper while moving towards the closest cyclotron 
resonance harmonic (cf. I) indicating considerable reduction 
of the phase. These observations confirm that at high intensi- 
ties microwave photoresistance can no longer be described by 
Eq.dB. 

Using the data in Fig.[T]we extract the amplitude A and the 
phase 6 for n — 2 and present the results in Fig.|2a) and 
Fig.|2jb), respectively, as a function of microwave power nor- 
malized to its maximum value, P/Pmax- Plotted in such a 
way, the data clearly demonstrate the existence of two dis- 
tinct regimes. At low intensities, the amplitude grows roughly 
linearly with power and the phase is power independent. 
Above a certain intensity, marked by I, the amplitude be- 
comes strongly sublinear and, simultaneously, the phase starts 
to decrease. 

The crossover between the regimes of low and high intensi- 
ties can also be established by examining the evolution of the 
oscillation waveform. Indeed, as illustrated in Fig.|2jc)-2(e) 
showing photoresistivity Spu:{e) at 20 dB, 10 dB, and dB, 
respectively, the waveform undergoes dramatic changes. At 
low intensity [Fig.|2jc)] the oscillations are well described by 
a damped sinusoidal, in agreement with Eq. ([T]i. However, al- 
ready at intermediate power [Fig.|2d)], one clearly sees that 
the extrema are pushed towards integer e. At maximum power 
[Fig.|3e)], the extrema become very sharp reflecting dramati- 
cally reduced phase. 

To explain the observed crossover we construct a theory 
based on quantum kinetics approach '^'^^•^-^ for arbitrary V. 
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FIG. 2: (Color online) (a) Amplitude A and (b) phase 5 measured at 
the second MIRO maximum versus normalized microwave intensity 
P/ Pmax- Solid lines in (a) represent linear and square-root depen- 
dences (see text). Photoresistivity 5pij{e) at attenuation of (c) 20 dB, 
(d) lOdB, and (e)OdB. 



Our main result is given by 

PD 

where, for circular polarization, 

J-(6 



2A2[-F(6)-1], 



267(07(0' 



-1 , -1 



7(6- 



(3) 



(4) 



^ = 2a/P sinTre, r^^ {t^) is flie transport (inelastic) scatter- 
ing rate, and the prime denotes the derivative with respect to 
e. 

For the model of mixed disorder^ with sharp (smooth) dis- 
order scattering rate r^y^ (Ts7^) and a typical scattering angle 
off smooth disorder x^/^ ^ 1, we find 

, 1 



7(0 
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'^sh 



1 - xO/2 



rs,n (1 + X^2)5/2 



(5a) 



.(5b) 



In the limit of smooth disorder, rj^ — 0, the first term in 
Eq. (in reproduces the displacement contribution obtained in 
Ref. 62 while the second term agrees with the inelastic con- 
tribution obtained in Ref.l63l assuming V <C t/tq. However, 
the real situation in the smooth disorder limit is believed to be 
more complicated. Already in the regime linear in microwave 
power, contributions from additional quadrupole and photo- 
voltaic mechanisms, due to the excitation of second angular 
and first temporal harmonics of the distribution functions, re- 
spectively, become relevant, and at higher power the num- 
ber of important angular and temporal harmonics increases.— 
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FIG. 3: (Color online) (a) Amplitudes Ai{P) and (b) phases Si{V) 
calculated using Eq. ^ with x ~ 0.022, to — 20 ps, r — 0.86 ns, 
and Tin = 0.10 ns. Asymptotes (lines) cross at V*^. [J-{e) — l]/27re 
vs e calculated for (c) P = 0.1, (d) P = 1, and (e) V = 10. 

In this work we consider an experimentally relevant situa- 
tion when the disorder contains both smooth and sharp com- 
ponents. The latter is assumed strong enough so that its contri- 
bution to the transport relaxation rate is comparable to that of 



the smooth disorder, t 



> 



In this case the quadrupole 



and photovoltaic mechanisms become negligible^^ while the 
displacement and the inelastic contributions can be separated 
[see Eq. (01)]. We analyze our results in terms of the reduced 
amphtude, 



A{V) 



27re+ 



(6) 



and introduce indicies i = sm, sh, in to describe smooth dis- 
placement, sharp displacement, and inelastic contributions, 
respectively. At low V, we find, in agreement with Eq. ^ 
and Refs.lSllillil 



(7) 



where 77^ 



Gx^T/Tsm, Vsh = 3t/2ts/i, and r]^n = 2tj„/t. 



In the limit of high V, all contributions show a square-root 
(inverse square-root) dependence of the amplitude (phase), 

Air) ^ v^(v:vy^^ s,{r)^{vt/v)'/\ (8) 

Here, Vt„, ~ 0.05/x, V^n ^ 0.4, Vt^ ~ min{l, r/2T„}.™ 
Calculated dependences of the amplitudes, Ai, and phases, 
6i, on microwave intensity V are shown in Fig.[3la) and 3(b), 
respectively. Distinct crossovers in both Ai and 6i depen- 
dences are evident for all three contributions. In Fig.[3lc)-3(e) 
we present [J-{e) — l]/27re, containing all contributions, as a 
function of e for (c) V = OA, (d) V = 1, and (e) 7^ = 10 
demonstrating considerable reduction of the phase at higher 



V, much the same as in the experiment, Fig.|2c)-2(e). While 
our data is in a good qualitative agreement with the theory, 
the observed crossover is noticeably sharper than theoretically 
predicted [cf. Fig.|3la) and Fig.|2a)]. Possible reasons for this 
discrepancy are heating, uncertain polarization, and the inho- 
mogeneity of microwave field, which are not included in our 
theory. 

We now illustrate the appearance of the crossover for the 
sharp displacement contribution given by the first terms in J^, 
Eq.® and in 7(^), Eq. (gbli. Since |C(e)| < 2VV, at low 
intensity we can use the Taylor expansion of 7(^) leading to 
Eq. ([T]l. At high intensity, the variable ^(e) spans an inter- 
val wide enough for the amplitude to be determined by the 
first maximum of the derivative dj{£^)/d^. As the first term 
in Eq. ( l5bb involves Bessel functions, its variation scale is of 
the order of unity. For that reason the maximum occurs at 
C+ = - 1 yielding e+ ~ n - l/27rVP. The amplitude 

then scales as (i^+/(ie+ w 2tt\/V, entering the full derivative 
rf7(^+)/de+. Explicit expression for the sharp disorder con- 
tribution is obtained substituting ( Bbl i with t^^^I = to the first 
term of Q, 



Tsh ^ 

+ 2^eV^cos(^e) Ji(0 [J2(0 " 3Jo(0] }• (9) 

Since the second line of Eq. (|9]) contains a factor 27re, it dom- 
inates at e > 1 and gives rise to the crossover as discussed 
above. Similar considerations apply to smooth displacement 
contribution, given by the second term of Eq. (Bbl l. 

The inelastic contribution, given by the second term of 
Eq. (|4|i, represents the balance between the excitation (numer- 
ator) processes, driving the distribution function out of equi- 
librium, and its relaxation (denominator). At Tin ^ the 
crossover occurs at V* ~ 1, very much like the sharp dis- 
placement contribution discussed above. Indeed, the excita- 
tion rate reaches its maximum at ^ = ^+ ~ 1 and, since 
7(^+) w I/tq, the relaxation rate remains roughly constant, 
sa T^^. In the opposite case of Tin ^ it is 7(^) in the relax- 
ation rate, varying in a narrow interval |^(e)| < \/t jTin <C 1, 
which controls the crossover. As a result, f ^ 1, and the 
theory of Ref. is applicable with the result 



47r7^esin27re 



■/Ti, 



2V sin^ 7re 



(10) 



The result (fToll holds provided V sin^ (7re) ^ t/tq — X^^ 
X <C 1 for the case of smooth disorder, t^^ — 0. However, 
when T~^ > XTrm the applicability range of Eq. ( fTOb is re- 
duced to V sin^(7re) ^ 1, as follows from Eqs. (|4|i and (|5]). 

Deviations from the linear power dependence of the numer- 
ator and the denominator are due to multiphoton processes and 
are captured by these equations. The extremum of Eq. ( fTOb is 
reached at V sin^(7re+) = r/2Ti„, and in the regime of slow 
inelastic relaxation Tin ^ the one-photon approximation 
suffices to describe the maximum of MIRO adjacent to the 
cyclotron resonance harmonics. Away from e ~ n and at high 



4 



intensities the inelastic contribution is strongly suppressed be- 
cause each impurity scattering event contributes to the relax- 
ation which saturates at I/tq, 7(^) w 0. This occurs because 
of the multiphoton processes which we discuss below. 

Observed crossovers reflect the transition from single- to 
multiphoton regime at elevated intensities. The rate of iV- 
photon scattering process is J^(-\/27'(l — cos6))/Tg, where 
JAf(a;) is the Bessel function of the order N and Tg'^ is the 
scattering rate on angle 9.— Since Jn{x) is strongly sup- 
pressed for N > X, a typical number of participating pho- 
tons is Ni w max{l,6i\/V}, where 6i is the contribution- 
dependent characteristic scattering angle. As a result, multi- 
photon processes become relevant at 7^ > V* w 1/^1- Since 
the scattering angles off smooth and sharp disorders are para- 
metrically different, so are the above crossovers. 

The square-root scaling of the amplitude, Ai cx Ni cx \/V, 
atV > V*, is naturally explained by noting that the fraction 
of electrons which can absorb or emit A'^^ -photons is given by 
Nifku/ep, where ep is the Fermi energy. The inverse square- 
root dependence of the phase St ~ {^Ni)^^ cx V^^l"^ follows 
directly from Nie^ = n — 1/4. 

In summary, our experimental and theoretical studies of mi- 
crowave photoresistance of a very high-mobility 2DEG over 
a wide range of microwave intensities revealed two distinct 



regimes. In one, low intensity regime, the oscillation am- 
plitude grows linearly with power while the phase remains 
power independent. In the opposite, high intensity regime, 
the amplitude exhibits square-root dependence of the ampli- 
tude and the inverse square-root dependence of the phase. The 
details of the crossover between these two regimes depends 
on the ratio Tm/r. For slow inelastic relaxation, Tm/r > 1 
the crossover is described by Eq. (fTOli.^^ For Tin/r < 1 the 
crossover occurs when a typical number of emitted (absorbed) 
photons becomes large, and is governed by multiphoton pro- 
cesses. Since the corresponding crossover scale depends on 
the type of disorder and the transport mechanism it can be 
potentially useful in separating different contributions. Taken 
together, our findings demonstrate an important role of multi- 
photon processes at elevated microwave power helping to un- 
derstand prior experiments-'^''^^"^" reporting sublinear power 
dependence of the oscillation amplitude. 
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The power dependence of the phase was investigated in Ref-pTOl 

reporting an intensity-independent value, 5 ~ 1/4. 

The inelastic mechanism has been investigated only in the regime 

V -C t/tq.^^ It has also been shown that in the case of smooth 

disorder different contributions cannot be separated and the 

MIRO undergoes several crossovers with increasing microwave 

intensity.'* 

The growth of the zero-field resistivity indicates an increase of 
electron temperature by about 1 K from 5 dB to dB. Due to such 
heating (and relatively low microwave frequency) zero-resistance 
states are not observed. 

These contributions determine the oscillatory part of the Hall con- 
ductivity for both types of the disorder.— 

Using X = 0.022, tq = 20 ps, r = 0.86 ns, and r,„ = 0.10 ns. 



we estimate P* 



1 P* 



2.4. 



In this case the temperature dependence of MIRO amplitude is 
expected to cross over from quadratic, oc to linear, oc T. 



